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I . INTRODUCTION 


Since the introduction of plastic filna in the nanufacture of high 
altitude halloona, acientleta have been able to exploit this type of plat- 
form in their studies of the atmosphere and space. The engineering commun- 
ity has attempted to keep pace with the demands of the experiments by pro- 
viding systems that will carry heavier loads to higher altitudes for longer 
periods of time. Due to the enormous value of these payloads, the relia- 
bility of balloon vehicles must be very high. Although many improvements 
have been made in the areas of film manufacturing, balloon fabrication and 
operations, there is a need for Improved design and analysis techniques. 

It is believed that the failure of balloons is directly related to the fail- 
ure of the film in service. Whether these failures are caused by overstress 
of good quality film or due to flaws or damaged film, the failure is driven 
by the state of stress In the film. This report is intended to address the 
nature of the stress state In a film rather than tie failure Itself. 

In order to improve the designs of scientific balloons, the ability 
to determine the state of stress (and strain) in the balloon throughout the 
entire flight regime must be developed. This necessitates Improved capa- 
bilities in three areas: the measurement of actual strains in balloons dur- 

ing real or simulated flights, the ability to calculate the stress field 
from the measured strains, and the analytical capability to predict the 
stress and strain in a balloon for all flight conditions. It is the long 
term aim of the current program at SwRI to satisfy these three requirements. 

The most urgent need of these three is the development of a device 
for accurately measuring stress or strain in balloon film since none current- 
ly exists. In a separate paper [1], Rand documents the development of a thin 
film strain transducer which promises to satisfy this need. The new gage has 
a low effective modulus so as to not Interfere with the strain field in the 
film, is highly sensitive to strains in the longitudinal direction and rela- 
tively insensitive to strains in the transverse direction, and is compensated 
for the extreme temperature variations common to balloon flights. 

With this new transducer, it should be possible to measure strains In 
the film at various locations along the surface of the balloon. In order to 
obtain stress values from these strain measurements, an appropriate consti- 
tutive equation must be developed which %d.ll take into account the interre- 
lationship of stress, strain, temperature and time. This material model 
would not only be usefxil in the analysis of flight data but also to improve 
existing balloon analysis techniques by incorporating a more realistic ma- 
terial model. 

The objective of the research documented in this report was to develop 
a constitutive relationship for polyethylene film which would allow the pre- 
diction of stress in the balloon film if the strain history and temperature 
history were known. It was determined that polyethylene could best be model- 
ed as a nonlinear viscoelastic material, which luplied that stress is a func- 
tion of the material's history. This realization led to a concerted effort 


to develop a numerical technique which would allow prediction of the stress 
in an efficient manne,*. To this end, an effort was made to combine pre- 
viously obtained test data with viscoelastic theory to develop a material 
model of the polyethylene film and a numerical <:echnique which can apply 
this model to the prediction of stress in the film at any time. The report- 
ed work concentrated on the case of uniaxial loading because of its simpli- 
city and availability of test data. F^rthar theoretical and experimental 
studies are required to develop a working model of the material for uniaxial 
loading and to generalize the results to the case of biaxial loading. 
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II. VISCOEIASTIC THEORY 


During the course of this program, It became apparent that polyeth- 
ylene film is a highly nonlinear viscoelastic material. The equations 
needed to describe the material respotse at a given temperature are complex 
with numerous functions which must be evalu'ited through testing. In this 
section, the general equations used to represent a nonlinear viscoelastic 
material are developed. Initially, a brief discussion is given of linear 
viscoelastic theory to acquaint the reader with this area. Data from actu- 
al tests on polyethylen: film are presented to document the existence and 
extent of material nonlinearity. Finally, a complex constitutive relation- 
ship is developed which will be used to model polyethylene film in the next 
section. 

A. Viscoelasticity 

Many solids, including those known as plastics, are referred to as 
being vlscoex^stic because they exhibit flow in addition to their elastic 
characteristics. This flow characteristic is often referred to as "creep" 
and is readily seen in a simple experiment. In this test, a length of 
material is "instantaneously" loaded in uniaxial tension with a constant 
force, as seen in Figure la. The resulting strain contains an initial 
elastic component plus a component which increases with time (creep;*) as 
seen in Figure lb. 

The creep compliance is defined as the time-varying ratio of the 
strain and stress as given by Equation (1) 


D(t) 


fit! 

0 

c 


( 1 ) 


It should be noted that for an elastic material, strain would be constant 
for this test, so that the compliance would be a constant. 

A dual response of a viscoelastic material is known as stress relax- 
ation, which refers to the response of the material under a constant strain. 
This response is readily seen in a simple experiment in which a constant 
strain is applied "Instantaneously" to a length of material. The load re- 
quired to hold the specimen at the constant strain level is seen to decrease 
with time from the initial load. The strain and resulting stress are seen 
in Figure 2. 

The relaxation modulus is defined as the time-varying ratio of the 
stress and strain as given in Equation (2) 

E(t) - 4^ (2) 

c 
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Per an elastic material, the stress would be constant, and Uie modulus 
would be comparable to Young's w>duluc, a constant. It should also be 
noted that the compliance, D, and the modulus, E, are reciprocals for an 
elastic material. 

Because of the wide range of times covered in a creep or relaxa- 
tion teat, the creep compliance and relaxation modulus are normally dis- 
played on log- log curves. Examplee of these properties are given In 
Figure 3. 

B. Linear Viscoelastic Theory 

From the definition of relaxation modulus. Equation (2), If a con- 
stant strain, Is applied at time - 0, the stress will be given by the 
relation 


o(t) - E(t) 


(3) 


For a linear viscoelastic material. If t%K> strains are superimposed and 
applied at time “ 0, the stress will simply be proportionately larger md 
given by the relation 


o(t) - [e^ + e^l E(t) 


(4) 


However, for the case In which is applied at t ■ 0 and is applied 
at time, t - t^, (Figure 4), the total stress at some time after tj^ may 
be approximated by the relationship 

o(t) ■ c E(t) + c E(t - t, ) (5) 

Tills theory of linear superposition of the two stresses at their elapsed 
times Is referred to as Boltzmann's superposition principle. For a gen- 
eralized strain history of the form shown in Figure 5, the stress may be 
approximated by the relationship 

N 

o(t) ■ ^ 

1-0 

For a continuously varying strain history. Equation (6) may be rewritten 
as 
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Strain Strain 



Tine Tine 


Figure 4. Superposli.lon of Sci esses 



Tlon Tine 


Figure 5. Generalized Stress-Strain History 
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( 7 ) 


t 

a(t) • J“ E(t - t) -|^ dx 

o 


It Is now apparent that for a linear viscoelastic material, thw stress at 
any tiue can be calculated approximately using Equation (7) If the strain 
history is known and a relationship exists for the relaxation modulus. 

By a similar derivation, a relationship for strain can be developed cor- 
responding to Equation (7) 


t 

e(t) » D(t - T) II dx (8) 

o 


The relaxation modulus and creep compliance for any material must 
be determined experimentally. The simplest methods are the creep test 
(constant stress) and the relaxation test (constant strain) described 
earlier. The major drawback to these tests is the inability to experi- 
mentally apply a finite stress or strain Instantaneously. This results 
in errors in nmasurement for short times. Thus, these tests are normal- 
ly used to obtain values of the modulus and compliance applicable for long 
times only (times greater than a minute) . For short time data, a more ac- 
curate value of the relaxation modulus of a linear viscoelastic material 
can be obtained from constant strain-rate tests. 

The constant strain-rate test consists of applying a constant rate 
of elongation and recording the load history. The strain history is taken 
to be 


e 


o t < 0 

< 

Rt t > 0 

V ~ 


(9) 


where R is the constant rate of strain. The stress history, o(t), is ob- 
tained from the loading data. Using this information and Equation (7) , 
one can develop a relationship for the relaxation modulus over a limited 
time span. Specifically, substituting Equation (9) into Equation (7) 
gives the result 


t 

o(t) = R f E(t - x) dx (10) 

o 
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By chttiging the variable of integration 




V ■ t - T, dv ■ -dx 
Equation (10) takes the form 

t 

o(t) - R J" E(v) dv (11) 

o 


This equation can be solved for B(t) by differentiating 


t 

o(t) « E(v) dv 
o 


1 . 

R dt 


E(t) 


But R dt > de, resulting in the relationship 


E(t) 


do 

de 



( 12 ) 


Thus, it is seen that the relaxation modulus at the time t - e/R is given 
by the tangent modulus (slope of the stress-strain curve). Using this 
technique, the relaxation modulus can be obtained over the time-span from 
fractions of seconds to minutes by conducting multiple constant rate tests 
at a variety of strain rates. 

Once material characterization has resulted in the derivation of the 
relaxation modulus, then the creep compliance can be derived from it, and 
vice versa. This can be seen by considering Equation (7) for the case of 
the strain history, e « D(t). For this case, the predicted stress is a 
unit step stress, such that Equation (7) becomes 
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1 


(13) 



E(t - t) 


3D 

77 


dt 


o 


Likewise, Equation (8) may be transformed to the form 


t 

1 - D(t - t) ~ dx (14) 

o 


Thus, once E(t) Is known, one can theoretically solve either equation In 
terms of D(t), ^nd vice versa. 

C. Relaxation Modulus of Polyethylene Film 

Constant strain rate data was generated for polyethylene film at a 
variety of rates and temperatures [2]. Curves for the relaxation modulus 
were generated from this data usli;g Equation (12) . Figure 6 displays the 
stress-strain curves and the derived relaxation modulus for polyethylene 
film tested at room temperature (T » 23*C) and three different rates of 
strain. The deviation of the moduli derived for each strain rate In Fig- 
ure 6b from a single curve was caused by -material nonlinearity. 

To this point all theoretical considerations have been based on 
the as;^jmption that the material behaves in a linear viscoelastic manner, 
l.e., that the relaxation modulus, E, and the creep compliance, D,-are 
functions of time (and temperature) only and not the level of strain. 

Since both the time dependence and the nonlinearity of the material pro- 
duce curvature in the constant strain rate data displayed in Figure 6a, 
the existence and extent of nonlinearity could not be determined by simple 
inspection of this figure. However, the extent of nonlinearity can be 
seen plainly in Figure 6b by observing the dashed lines which were obtain- 
ed by connecting the moduli measured at the same strain levels. The modu- 
lus is seen to decrease dramatically with increasing strain. 

For a nonlinear viscoelastic material, Boltzmann's superposition 
principle, as represented by Equation (5), is no longer valid. For this 
reason, our single integral representations of the modulus. Equation (7), 
and the compliance. Equation (8), are no longer valid. 

D. Nonlinear Viscoelastic Theory 

Various theories have been developed for representing the constitu- 
tive relationships of nonlinear viscoelastic materials [3]. The simplest 
and most tractable of these appears to be the theory of Schapery [4, 5] 
which is a modification of the single integral solution represented by 
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Equations (7) and (8) for constant teaperature. This theory takes Into 
account the obaervatlou that the nodull at hi^ strains in Figure 6b 
(dashed lines) are the same as those at low values, except for horizontal 
and vertical translations. Mathematically, Schapery rewrites Equations 
(7) and (8) for linear viscoelastic response 


t 

o(t) - E^ c(t) + f AE(t - t) dx (15) 

o 


t 

e(t) - o(t) + f AD(t - T) l^dx (16) 

o 


where £» represents the equilibrium or final (long time) value of the 
modulus. Do represents the initial value of the compliance, and 


AE(t) - E(t) - E 
00 

(17) 

6D(t) - D(t) - D 

o 

(18) 


Schapery extended these equations to nonlinear response on the basis of 
his theory of the thermodynamics of irreversible processes. By assuming 
certain simple forms for free energy and entropy production, he derived 
the constitutive equations for constant temperature 


o(t) 


h 

00 


E e(t) + h, 

00 1 



where 


P 


t 

p(t) - J' 

o 



(a^ > 0) 


(19) 


( 20 ) 
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T 

p' - p'(t) - r ^ (21) 

i ' 

- a^Ie(t')] 
and 


c(t) - gQ o(t) + Si f - 


3 8-0 

r) dx 


wliere 


- r ^ > 0) 

y a o 


r - *'( t ) - / ^ 


- a^[o(t')] 
o o 


(22) 


(23) 


(24) 


In these relationships, hao and h]^ are functions of c(t), h2 and a^ 
are functions of e(x), go and are functions of o(t), and g 2 and a^ are 
functions of o(x). For Equation (19), ho,, h^, ho, and a^ tend toward unity 
as the strain approaches zero, so that Equation (19) reduces to the linear 
relationship. Equation (7). For Equation (22), gQ, gj, g 2 , and a^ tend 
toward unity as the stress approaches zero, so that Equation (22) reduces 
to the linear relationship. Equation (8). 

Bepresentatlons of the nonlinear modulus and compliance are obtained 
by substitution of a constant strain, e. Into Equation (19) and a constant 
stress, o. Into Equation (22) «ihlch yields 
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B 

n 



( 25 ) 


(26) 


where Bq and Dq are Che nonlinear modulus and nonlinear compliance respec- 
tively. 

E. Temperature Dependence 

All discussion to this point has assumed that the temperature remains 
constant. With an Increase In temperature the creep and stress-relaxation 
processes are accelerated, resulting In a large variation of the relaxation 
modulus and creep compliance with temperature for most polymeric materials. 

It has been obeerved for these materials that the shapes of the modulus- 
time and compliance-time curves (Figure 3) do not change radically, but do 
suffer a translation parallel to the log t axis. 

For a significant number of materials, this shift of property functions 
parallel to the log t axis without change Is sufficiently accurate to be the 
basis of a thermo-vlscoelastlc theory. A material which exhibits this behav- 
ior Is often termed thermorheologlcally simple and the behavior Is termed the 
time-temperature equivalence. An example of time-temperature equivalence Is 
given in Figure 7. 



Figure 7. Time-Temperature Dependence 


Norland and Lee [6] have shown that for a linearly viscoelastic, 
thermorheologlcally simple material, the constitutive relationship of Equa- 
tion (7) takes the form 


14 


(27) 


t 

o(t) - / E, (C - t') ^ dx 
o 

where 





reduced tine 


(28) 





(29) 


and 


- a {T(t)} 


In Equation (27), EB(t) la the stress relaxation function at a specific 
reference tenperature, Tg|, which is temed the "master" curve; a^ is the 
shift function such that 


*T 


T - T 


1.0 


For the case of the temperature held constant at X ^ 1^, Equations 
(28) and (29) can be simplified 




*T 


(30) 




“X 


(31) 


and Equation (27) takes the form 


15 


0(t) 


(32) 



It should be noted that 


E (t. T) - E_ (O 


(33) 


That is, the relaxation modulus applicable for the temperature T can be 
obtained from the master curve by shifting the time axis according to 
Equation (30) or 


log C * log t - log a^ 


(34) 


iihere log a^ corresponds to the shift along the time axis in Figure 7. By 
convention 


a 


T 


> 1 for T < T 
m 

< ■ 1 for T ■ T 
m 

< 1 for T > T 
V. m 


(35) 


When the temperature field varies with time, the reduced times, C and C', 
are given by Equations (28) and (29) which take into account the entire 
past history of the temperature field. 

A very Important feature of the time-temperature equivalence is the 
ability to use data obtained at different temperatures to derive the relax- 
ation modulus at any given temperature. This is demonstrated in Figure 8 
in which relaxation moduli obtained at various temperatures using constant 
rate tests were used to construct a master curve extctndlng over an extreme- 
ly large time span. 

To this point, the discussion on temperature dependence has assumed 
that the material is linearly viscoelastic. This is necessary for Equations 
(27) and (32) to be applicable. However, the time-temperature equivalence 
has been seen to hold for nonlinear viscoelastic materials as well. Equa- 
tions (33) and (34) can be applied equally mil to the nonlinear equations 
derived in the last section. In this case. Equations (19) and (22) may be 
rewritten 
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data master curve 



Figure 8. Construction of a Stress Relaxation Master Curve 
of an Anorphous Polymer (Reference 7) 


t 

Eg-h.E. + hi / me 

O 


3h,e 

6-) d. 


where 


t 



a 


cT 


a^^ (c(t), Tj 


(36) 


(37) 


(38) 


and 


c(t) 


80 Do °<t) + 


/ 


AD(y - Y') -37- dx 


(39) 


17 


whare 


Y 


Y 



t 



(AO) 


(Al) 


•oT " ‘oT 


The najor change In thaaa equatlona waa the addition of tanperature depend* 
anca to a^ and ao defined earlier. 

Finally, It haa bean noted [8] that the condition of tiBe-tanpara* 
tura equivalence is not a univaraally valid rule. In order to Increaaa its 
range of applicability, atte^its have bean nada to nodal data by Including 
a nagnification factor (a vertical riiift of the curvea parallel to the loga- 
rithadc nodulus axis) In addition to the horizontal shift already discussed. 

Such horizontal and vertical shifts have been used previously [8] to 
develop neater curves for several types of polyethylene filn. The naster 
curve for Hlzex, displayed in Figure 9b, was fomed by shifting the stress 
relaxation curves of Figure 9a. The shift factors were 


horizontal log a^ ■ *0.035 (T * 20*C) 

vertical log h^ - -0.0083 (T * 20*C) 


The presence of the vertical shift iaplies that this polyethylene filn (Hlzex) 
is not a themorheologlcally-slnple naterial. However, the techniques pre- 
sented above appeared to be sufficient for developing a realistic naster curve 
for this naterial. 
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(a) 



(b) Master Curve 


Figure 9. Streaoi Relaxation Curves of Polyethylene (Hizex) 
at Various Teaq>eratures [8] 


19 




III. SIMPLIFIED MATERIAL MODELS 


Tha purposa of this ongotog program Is to davalop tha ralstlonshlps 
for charactarislng tha raspoasa of balloon film undar load. Of initial 
Importanca is a ralationship which would allow thr calculation of strass 
for a known atrain hietory. This is naadad to procass data attalnad using 
strain gagas during balloon matarial tasts and actual flights. 

Bacausa of the complexities involved in developing a nonlinear modal 
of film vhi :h accurately predicts response over a wide range of tamparatura 
and loadLsg conditions, emphasis has bean placed on developing slmpllfisd 
models Initially. In this section, a review of previous work will be given. 
An explanation of a simplified analysis developed In this program will also 
be given. This simplified model was developed to lend Insights into the 
nonlinear processes and to demonstrate the feaaibllity of properly modeling 
film response with a nonlinear model. Emphasis in present %fork has been on 
model development for uniaxial loading while biaxial effects will be con- 
sidered in later work. 

A. Background 

The mechanical response of polyethylene has been studied under a 
variety of loading and temperature conditions by a number of investigators. 
Host notable has been the data collected by Alexander and Walssman (9] at 
Stevena Institute of Technology. Okamoto [10). and Webb [2] at Texas A&M 
University. Unfortunately, raw data cannot be effectively utilized in any 
logical analysis procedure. Alexander characterized the film as a linearly- 
elastic, orthotropic material with material properties evaluated as a func- 
tion of only temperature. However. Vebb [2] experimentally showed that the 
mechanical properties are a function of time as well as temperature. This 
suggested that the material is viscous in nature, although not necessarily 
viscoelastic. 

It was suggested by Schapery [11] that if the data are restricted 
to very smell strains, then linear viscoelastic characterization might be 
appropriate. Thus. Webb [2] conducted an extensive effort to characterize 
polyethylene film as a linear viscoelastic material. Several simplifying 
techniques were incorporated into his analyses including a linear time- 
temperature shift, a simple power law representation of the tensile modulus 
and the assumption that the modulus and compliance were reciprocals. In 
this program, initial work concentrated on developing a relationship for 
the creep compliance which was a rur-:tion of both time and teqperature. 

The compliance was chosen since it can be developed from creep test data, 
which is relatively easy to generate and docs not require extensive test 
equipment. The relaxation modulus was then assumed to be the reciprocal 
of the derived creep compliance. The resulting relatitmship for the modu- 
lus was 
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E(t) 


(42) 


• 20,000 



whtn 


log - -0.U6 T ••• 2.668 (43) 


«ttd 


E<t) is in psl 
t is in Binutss 
T is In *C. 


This rslstionship was basad on tasts of thraa thicknsssas of balloon film 
(0.5, 0.7, and 1.0 sdl) with a trada nane of StratofUa. 

As a chock of this relationship, constant strain rata tasts ware 
conducted and the results were checked against predictions snde using 
Bqustion (42) and the convolution Integral, Equation (7). Tor this case, 
tte convolution Integral took the fora 


o(t) • R J 20,000 

o ^ 

Constant rate data mbs generated for three teaperatures Mid three rates. 
Hesults of the predicted and aeasured stresses are given in Figure 10. 

A poor correlation was found between the aeasured and predicted 
stresses, especially at higher strains. The reasons for this are rather 
strait t forward: the siaple power law representation o, the asterlel was 

Inadequate, asterlel nonlinearities of the asterlal were not considered, 
teaperature shifting was overly slaplified, and the assuaptlon that the 
coapllance and aodulus are reciprocals was invalid. 

B. Slaplified Analysis 

In order to properly nodal the rer«ponse of polyethylene film, an 
expression for the relaxation aodulus was developed which took the form of 
Equation (25) 
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Stress (HPa) 



E_ 


11. E. + hj hj 4E 

G 


with four functions, h„, hj^, h 2 , and a^, vhich oust be evaluated axperlnen- 
tally for any given tenperature. An axtenalve teat prograa Is required to 
accurately evaluate these functions over the range of temperatures of inter- 
est. Since data from this type of experimental program does not currently 
exist, an effort %ias mede to use existing data to develop a simplified ex- 
pression for the nonlinear modulus. Such an expression, though Incapable of 
modeling the response of the material under all conditions, mould he useful 
in evaluating the effects of material nonlinearities on the prediction of 
stress and strain. 

One approach to developing a sinplifled expression for the modulus 
is to modify Equation (7) such that 


t 

a(t) « E^ (t - t) dT (44) 

o 


where the linear modulus has been replaced by its nonlinear counterpart, 
which has the form 


■ "e ' (^) 


(45) 


This simplified relationship is similar to Equation (25) for the case of 
h» ■ h]^ h£ ■ he and a^ “ ag* Equation (45) can be rewritten 


log E^(t) 


log h^ 


-loEE(^) 


(46) 


If experimentally obtained modulus data Is plotted against t on double-loga- 
rithmic paper, curves at different strain levels can be superposed by trans- 
lation along the two axes. The amount of vertical shift corresponds to log 
hg and the amount of horizontal shift to leg Se. This process Is similar to 
that explained In Section II for temperature shifting. 

The only data currently available for evaluating the two nonlinear 
coefficients In Equation (45) is the constant strain rate data obtained by 
Webb [2]. The stress-strain data given In Figure 11 is for 1.0 mil Strato- 
fllm which was tested in the transverse direction. Information concerning 
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Figure 11. Uniaxial Stress-Strain Behavior of 1.0 Hil Stratofllm 




the relaxation nodulua can be obtained by measuring the slope of the stress- 
strain curve. According to Equation (12) » the relaxation modulus at the 
time t " e/R Is given by the measured tangent modulus (slope of the stress- 
strain curve). The measured values of the tangent modulus obtained from the 
data In Figure 11 for T " -23”C Is given In Figure 12. A separate curve was 
generated for each strain rate. 

It should be noted that these curves were generated from data for 
which the strain Is constantly changing, whereas the relaxation modulus Is 
defined as the ratio of stress to strain for a constant strain. Thus, to 
obtain curves for the relaxation modulus from tangent modulus curves, points 
on each curve in Figure 12 for a given strain must be connected. As seen in 
the figure, dashed lines have been added which connect the points of constant 
strain. These dashed lines represent the curves for relaxation modulus at 
strains of 1, 2, 3, and 4 percent. The fact that these curves are not super- 
posed indicates that the modulus Is a function of strain, and therefore, non- 
linear. As the strain level decreases, each of the tangent modulus curves 
tend toward a single curve, which represents the apparent linear relaxation 
modulus. 


Looking again at Equation (46), the linear modulus Is represented by 
E(t) and the dashed lines can be superposed on the linear modulus curve by 
vertical and horizontal translations, with the vertical translation repre- 
sented by log h^ and the horizontal translation by log ae. In order to sim- 
plify the model being developed, a decision was made that only vertical 
translations would be used to shift the data, that Is, ae » 1. A plot of 
the measured vertical translations are given In Figure 13. The expression 
for he derived from this curve took the form 


d log h^ 
de 


-17.2 


or 


h 


-n.lt 

e 


(47) 


This relationship, along with Equation (45), Is sufficient to model the ma- 
terial If the temperature remains at -23*C. However, a relationship is 
needed which can be used to describe the response of film to the entire range 
of temperatures experienced in balloon flight. Such a relationship can be 
developed by deriving relationships similar to Equation (47) for a range of 
tenq>eratures and then developing a master curve as described in Section II. 
The horizontal and vertical shifting done In constructing the master curve 
will result in Equation (45) taking the form 
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"• ^ (t) 


(48) 


E„(t) 


where It has been assumed that a^ ■ 1 for all temperatures. 

Expressions similar to Equation (47) were derived for much of the 
data generated by Webb [2], and is displayed in Figure 14. From this fig- 
ure, an expression for as a function of temperature was found 


^ (-27.63 + 0.394 T) c 

h ■ e 


(49) 


with T in *C. 

Curves similar to the linear modulus of Figure 12 were developed for 
several series of tests, including O.S mil film tested in the machine direc- 
tion and 1.0 mil film tested in tlte machine and transverse directions. 

These apparent linear moduli are plotted in Figure 15 for the four test 
temperatures. In a process similar to that demonstrated in Figure 8, a 
master curve wr.s 'constructed from these curves, and is presented in Figure 
16. Thus, the llnenr relaxation moduli obtained at various temperatures 
over a short time sp«.'n were used to construct a master curve extending over 
an extremely large time span. A modified power law was developed which ad- 
equately represents the master curve and which has the form 


E(t) 


0.0069 GPa 



2.30 GPa 


5 0.094 


(50) 


with t in seconds. 

The amount of horizontal and vertical shifting which was required to 
obtain Figure 16 is plotted in Figure 17. From this figure, the expressions 
for aj and hj were found to be appr :>ximately 


(1.635 - 0.0716 T) 

St - e 


(51) 


^ (0.409 - 0.0179 T) 

h^- e 


(52) 
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Figure 15. Derived Linear Relaxation Moduli for Stratofllm 
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Now, th« •xprasalona glvaa In Bquationa (48) through (52) ara aufflclant 
to daacrlba tha raaponaa of tha film ovar tha antira range of tanparaturae 
and tlnae which ara axpariancad in balloon fli^ta. 

Aa a partial chack of tha validity of tha alapla nodal juat devel- 
oped, a dacialon uaa nade to predict tha atraaa-atraln curvaa which had 
been uaad Initially aa Input data In tha derivation of tha a»dal. Thla waa 
dona by aubatltutlng tha alapllfled nonlinear nodal repraaantad by Equa- 
tlona (48) through (52) Into tha ralatlonahlp for atraaa as a function of 
atraln. Equation (44). Equation (44) waa nunarlcally Integrated for twelve 
caaaa coaprlaad of taata at four ta^»araturea at each of three atraln rataa. 
Tha raaulta of theaa calculatlona ara given In Figure 18. Conparlaona of 
thaaa raaulta with aona of the Input data given In Figure 11 ahows that the 
nonlinear nodal doaa a good job qualitatively of duplicating the atraaa- 
atraln curvaa over tha range of rataa and taeperatiuraa Involved. Thaaa pre- 
dict Iona greatly contraat thoaa nade ualng tha linear nodal developed by 
Webb [2] and dlaplayad ;>n Figure 10. The al^>llfled nonlinear nodal closely 
nodela tha decrease in nodulus with Increasing strain which la alnost con- 
plataly absent with the linear nodal. 

Although the slnpla nonlinear nodal shows good qualitative agreanent 
with the United data. It is still Inadequate as a nodal for quantitative 
calculations, as typified by Its inability to duplicate the naasurad re- 
sponse at T • ■f38*C. Farther %rork with this oodal has danonstratad its In- 
ability to properly nodal the case of cyclic loading, such as Is present in 
an actual balloon flight. 

In a slnpla cyclic loading test, a thin tensile film spaclnen vta 
subjected to the following strain history: 


■f2 an/nm/mln 
-2 nm/nsi/min 


i • < 


•f2 mn/ani/inln 


I -2 un/on/adn 
^■f2 na/nm/aln 


0 ^ t < 5 nin 
5 ailn £ t < 7 min 

7 nln £ t < 12 nln 
12 nln £ t < 14 nin 
14 min < t < 18 laln 


The tenperature was held constant at T ■ -f23*C. Typical results of this 
test are given in Figure 19. The predicted response of the film using the 
simple nonlinear model is given In Figure 20. An exanlnatlon of the pre- 
dictions show that the simple nonlinear nodel does not properly model the 
release process. This Is probably caused by tha inability to nodel all 
aspects of the material's history dependence with a single nonlinear func- 
tion, he* 
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Figure 18. Predicted Stress-Strain Curves for Stratoflla 
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Flgpra 19. Haasured Raaponae of ScratofllB to Cyclic Strain Hiatory 
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Figure 20, Predicted Reaponae of StratofUa to Cyclic Strain History 
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In the next section, four nonlinear functions are included in the 
model for the relaxation m^ulus. These additional terms should allow 
the model to more properly duplicate the actual response of balloon film 
to the fll^t environment. 
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IV. PROPOSED NONLINEAR MODEL 


In order to properly model polyethylene film, a more complex non- 
linear model must be used than the simple model considered in Section III. 
The desired model chosen for this study was based on the Schapery thermo- 
dynamic theory [4, 5]. The constitutive equation for the stress in terms 
of the strain history was given in Equation (19) 


r ah^c 

o(t) • h^ e(t) + J AE (p - p") di 

o 


with a corresponding nonlinear relaxation modulus given by Equation (25) 

E (t) - - h E + h- h- AE 

n e o# • 12 ^a^ j 


The functions h«>, hi, h 2 , and a^ are dependent upon the strain. More spe- 
cifically, hoD and hj^ are seen to be functions of the current strain, e(t), 
and h 2 and a^ are functions of the entire strain history, c(t). It should 
be noted that for a linear material, tv» > hi > h 2 ■■ a> <■ 1, and Equation 
(25) reduces to the linear relationship. Equation (17) 


E(t) - E + AE(t) 
00 


In order to evaluate Equation (25) for a given material, an extensive 
test program is required. In this process, relaxation tests are conducted 
at a variety of temperatures. For each temperature, a value is obtained for 
Ew, an expression is developed for the transient linear modulus, E(t), as a 
function of time, and expressions are generated for h«o, hi, h 2 , and ae as 
functions of strain. 

In order to generalize Equation (25) for all temperatures of interest, 
expressions are developed for the four nonlinear terms as functions of tem- 
perature as well as strain. This process is similar to that used to evalu- 
ate the nonlinear function, h^, of the simplified model as seen in Figure 13, 
with resulting relationships similar in form to Equation (49). Generalization 
of Ego and AE(t) in terms of temperature is accomplished in the manner dis- 
cussed in Section II. The linear modulus, Eq, -f AE(t), is plotted for each 
different temperature, and then shifted vertically and horizontally to form 
a master curve. The vertical shift is used to develop an expression, ’«■] , and 
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the horlsontel shift Is used to develop aj. The resulting expression for 
the linear portion of the modulus takes the form 


E(t. T) - + h^ AE (53) 


The full expression for the nonlinear modulus then has the form 


E„(t. T).|.^h,E. + hjhjhj4E(j-E-) (54) 

\ 6 T ' 


By cosi>lnlng several terms. Equation (54) may take the form of Equation 
(25) 


E^,(t. T) - h.' E. + hj 4E (^) (55) 


where h„ • h^ h «09 h^ • h^ h^t and a^ * a^ 

The key to developing this relationship for the nonlinear modulus Is 
the evaluation of the four functions hn, hi, h 2 , and a^ at each of the tem- 
peratures since the remainder of the process Is similar to that used in de- 
veloping the slnple model of Section III. These functions can be evaluated 
empirically through the use of a carefully conducted series of relaxation 
tests at a variety of temperatures. This process will not be explained. 

Initially, a series of one-step relaxation tests should be conducted 
at a variety of strain levels and temperatures. The strain histories will 
be of the form 


e 



t < 0 


0 < t < t 


For these tests, the proposed nonlinear relaxation modulus has the form 


E (t) 
n 


• “ h E + h, h- 

a eo 00 T. 2 
c 



(56) 
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for all tests at a given temperature. Now, the transient component of the 
modulus is seen to be 


■ “ (f ) 


(57) 


This relation may also be written 


log AE^(t) 


log hj^ h 2 + log AE 



(58) 


If the results of several relaxation tests w;;iiQucted at different levels 
of strain but the same temperature are plotted on a double logarithmic 
graph, they should appear similar to the dashed lines of Figure 12. Ideal- 
ly, one of the tests should have been conducted at sufficiently low strain 
to assure linear response. The results of the linear test would then cor- 
respond to the linear curve in Figure 12. If the relaxation tests were 
conducted for a sufficiently long time, values for the equilibrium moduli 
should be found from these curves. For the linear curve, this value cor- 
responds to Eos, and for the other curves, it will be heo E». Thus, rela- 
tionships can be obtained for both £<» and ho> for each temperature tested. 

A check of Eo, can be made by conducting a long term creep test with suffi- 
ciently low stress to insure linear response. The measured equilibrium 
creep value, D<», should be the reciprocal of Ea,. 

By subtracting the equilibrium values of the moduli, curves for the 
transient portions of the moduli, AEg(t), can be obtained. When these are 
plotted against time, t, on double logarithmic paper, it should be possible 
to superpose all of the curves by translating them along the two axes. 

This process was followed in the last section to obtain values for h^ and 
ae. In the present case, the horizontal shift will be equal to log and 
the vertical shift to log hi h 2 * Thus, if this process is successfully 
followed, values can be obtained for a^ and the product h] h£ for a given 
temperature. 

In order to obtain separate values for hj^ and h 2 , a series of two- 
step relaxation tests must be conducted at various ten^eratures. For these 
tests, the strain histories will be 
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t < 0 


(o 
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< c 
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b 


V 


0 < t < t. 

‘a i " S> 


Substituting this strain history into Equation (19) results in the follow- 
ing relationships [12] 


0 < t < t 
— a 


a ■ c 


a 


[' 


a 



(59) 


and 



(60) 


( 61 ) 
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The integrations used to derive Equations (59) and (60) were carried out 
recognizing that 3h2e/3t is zero except at t ■ 0 and t > t^. 

If one-step relaxation tests were conducted for the strains and 
e^, then values should already be known for Eg»f h«i^> has^t 
the products h^^ h 2 ^ and hi^ ho^. Thus, all terms are known in Equation 
(61) except for the ratio (hj^^/hj^^) which can then be obtained by matching 
Equation (61) to the two-step data at a convenient point (t > t^). At 
least one test must be mn in which was chosen to be sufficiently small, 
so that the material behaves linearly and h]^^ ■ 1. This allows evaluation 
of hj^A, and consequently, h 2 ^. In subsequent two-step tests, strains will 
be chosen for which have already been evaluated in terms of h]^ and h 2 , 
so that the value hj^^ will be known. Through the above process, values can 
be obtained for hoo, hi, h2, and a^ at each of the temperatures of Interest. 

A process has been described which should be sufficient to allow for 
the development of a constitutive relationship of the form of Equation (19). 
Tk>e tests required to evaluate this relationship completely are given in 
Tabic 1. This series of tests must be conducted at a sufficient nuad>er of 
temperatures, probably three or four, to allow for adequate evaluation of 
the constitutive relationship over the full range of temperatures Inherent 
in balloon flight. 


Table 1. Proposed Test Matrix For Each Temperature 


Test 

Nunber 

One-Step Relaxation Test 

Two-Step Relaxation Test 

a 

e 

a 

e 

b 

e 

1 




2 

^2 

"2 

"l 

3 


"3 

^2 

4 


^4 

"3 

5 

^5 

"5 

"4 


must be sufficiently small to allow only linear response. 
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V. DISCUSSION OF RESULTS 


The viscoelastic oature of thin polyethylene film has been care- 
fully detailed In the preceding sections of this report. Available data 
has been used to graphically demonstrate the nonlinear character of this 
material idilch has previously been modeled as either linearly elastic or 
viscoelastic. In addition, a test procedure has been described which will 
permit the evaluation of the various functions needed to produce an accu- 
rate nonlinear characterization. 

In order to assess the Importance of an accurate material charac- 
terization In the analysis of balloon stress and strain, a simplified com- 
parison of the results obtained to date and previous results Is revealing. 
Consider the problem of creep under constant stress conditions at -23°C. 

This vould be somewhat analogous to the circumferential stress state of a 
balloon under float conditions. Since laboratory data Is obtained In a 
relatively short period of time, then the corresponding material proper- 
ties such as those reported by Alexander [9] would predict relatively small 
strains. If a stress of 6.89 MPa (1300 psl) Is assumed then strain on the 
order of 1.5 percent would be predicted. However, a typical balloon must 
remain In this stressed state from 4 to 24 hours which will cause the strain 
to Increase to as much as 2.8 percent, an Increase of 86 percent from the 
original value. This estimate is based on the simplest linear viscoelastic 
model. Equation (42), which has already been shown to neglect the nonllnear- 
Itles due to strain. 

Another representative problem Is that of relaxation under constant 
strain. This would be somewhat analogous to the meridional stress state 
of a balloon since the load tapes will tend to restrain deformations in this 
direction. If a strain of 1.5 percent is assximed, then a stress of 6.89 MPa 
(1000 psl) would be predicted from short term data. This stress would relax 
to 3.62 MPa (525 psl) after 24 hours which represents a decrease of 47 per- 
cent from the original value. In a balloon, this would cause a load redis- 
tribution and the tape force would have to Increase to maintain meridional 
equilibrium. 

Another problem of Importance, especially during ascent. Is the tran- 
sient temperature history of the film. A method must be developed to take 
this event Into account. This may account for the observations by several 
Investigators that sunmer operations tend to have more apparent balloon 
failures than engineering flights conducted In the winter months. This may 
be caused by the substantially different strain distributions which would 
occur If a balloon Is held In the launch configuration at different tempera- 
tures prior to ascent. 

Each of these problems Indicate that transient response of the film 
must be understood If the state cf stress or strain is to be determined with 
sufficient accuracy to Identify causes of balloon failure. In order to ac- 
complish this, a selective set of laboratory tests directed by those respon- 
sible for the characterization Is necessary to obtain realistic values of 
stress from strains measured In flight. 
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